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Abstract. Significant progress in the description of quasar variability has been recently made 
by employing SDSS and POSS data. Common to most studies is a fundamental assumption that 
photometric observations at two epochs for a large number of quasars will reveal the same statistical 
properties as well-sampled light curves for individual objects. We critically test this assumption 
using light curves for a sample of ~2,600 spectroscopically confirmed quasars observed about 50 
times on average over 8 years by the SDSS stripe 82 survey. We find that the dependence of the 
mean structure function computed for individual quasars on luminosity, rest-frame wavelength and 
time is qualitatively and quantitatively similar to the behavior of the structure function derived 
from two-epoch observations of a much larger sample. We also reproduce the result that the 
variability properties of radio and X-ray selected subsamples are different. However, the scatter 
of the variability structure function for fixed values of luminosity, rest-frame wavelength and time 
is similar to the scatter induced by the variance of these quantities in the analyzed sample. Hence, 
our results suggest that, although the statistical properties of quasar variability inferred using two- 
epoch data capture some underlying physics, there is significant additional information that can be 
extracted from well-sampled light curves for individual objects. 
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INTRODUCTION 

Significant progress in the description of quasar variability has been recently made by 
employing SDSS data. Vanden Berk et al. (Ill], hereafter VB) compared imaging and 
spectrophotometric magnitudes to investigate the correlations of variability with rest- 
frame time lag (up to 2 years), luminosity, rest-frame wavelength, redshift, the presence 
of radio and X-ray emission, and the presence of broad absorption line outflows. Vari- 
ability on longer time scales was studied by de Vries, Becker, White, and Loomis 
hereafter dVBWL; yfl) who compared SDSS and POSS photometric measurements. 
Ivezic et al. ((4fl, hereafter, 104) used repeated SDSS imaging scans, which increased the 
measurement accuracy for magnitude differences by about a factor of 3-4 compared to 
studies based on spectrophotometric magnitudes, and also enabled analysis of the u and 
z band variability. 

Common to these and similar studies is a fundamental assumption that photometric 
observations at two epochs for a large number of quasars will reveal the same statistical 
properties as well-sampled light curves for individual objects. We critically test this as- 
sumption using light curves for a sample of ~2,600 spectroscopically confirmed quasars 
observed by the multi-epoch SDSS stripe 82 survey. We describe these observations and 
our methodology in the next section, and analyze the data in Section 3. 




FIGURE 1. (Figure 1 from 104) The left panel displays the measured structure function as a function 
of rest-frame time lag and wavelength (each pixel contains ~ 25 objects). The middle panel shows the 
structure function in a narrow range of rest-frame time lag and wavelength (100 days< AtRp <150 days, 
1000A< Xrf <2000 A), as a function of redshift and luminosity. The lines of constant variability are 
nearly parallel to the redshift axis, suggesting that the dependence of variability on luminosity is much 
stronger than the dependence on redshift. The right panel displays the difference between the data shown 
in the left panel and the best-fit model described by Eq. (Q]). The over-plotted line shows the composite 
quasar spectrum from [7]. The negative residuals at ~2800 A are due to the anti-correlation of continuum 
and Mgll line variability. 



SDSS STRIPE 82 DATA AND METHODOLOGY 

We study a sample of about 2,600 spectroscopically confirmed quasars (see [5] for the 
SDSS Quasar Catalog) imaged about 50 times on the average over 8 years (for details see 
[@]). These data were obtained in yearly "seasons" about 2-3 months long. We average 
observations within each season, and construct a structure function (see dVBWL for 
definition) for each object for an observed time lag of 1 year, and in each of the five SDSS 
bandpasses (ugriz). The total number of structure function data points for all quasars and 
bands is 10,370. We correct for cosmological time dilation and length contraction and 
obtain a fairly well- sampled plane of the rest-frame quantities A^ and Xrf (spanning 
from about 100 to 300 days and 1000 to 6000 A). 

To compare the statistical properties of quasar variability inferred from two-epoch 
data, and those based on light curves for individual objects, we use a quasar variability 
model from 104. Using ~66,000 magnitude difference measurements, Am, for a sample 
of ~13,000 quasars, they found that Am follow an exponential distribution, p(Am) 
exp{— |Am|/A c ), where the characteristic variability scale, A c is a function of rest-frame 
time lag (AtRf, days), wavelength {Xrf, A), and absolute magnitude in the i band 
(Mi) (the variability scale is related to the more commonly used structure function by 
SF = \^2A C ). Their simple model for the variability scale, 

SF= (1.00 ±0.03) [1 + (0.024 ±0.04) M t ] (1) 

describes Am measurements to within the measurement noise (~0.02 mag). Note that 
there is no dependence on redshift (see Fig.Q]). 



We attempted to fit the same functional form to the SF data obtained for individual ob- 
jects from stripe 82. We find that the observed SF is still well-fit using a single exponent 
for Xrf and AIrf, except with a value for the exponent of 0.47 ± 0.02. The parame- 
ters used to model the observed SF are measured in the ranges 1000 A < Xrf< 6000A, 
100 days < AtRF < 300 days (compared to 0-600 days for the two-epoch sample), and 
—29 < Mi < —23. This modified model SF provides a good description of the observed 
SF; the ratio of the observed to model SF does not show any systematic behavior with 
respect to At^p, Xrf, Mi, or redshift. 

The similarity of the mean structure function computed for individual quasars and the 
structure function derived from the two-epoch observations of a much larger sample 
suggests that the statistical properties of quasar variability inferred using only two 
epochs for a large number of quasars are representative of the underlying physics. This 
is at least true for the overlap in their timescales (i.e., for timescales ~300 days or less). 
We proceed with analysis of the distribution of the ratio of observed (SF £, S ) and modeled 
{SF mo dei) structure functions. 

ANALYSIS: THE STRUCTURE FUNCTION VARIANCE 

The distribution of SF t, s /SF m0C { e i values for the two-epoch sample analyzed by 104 is 
shown in the upper left panel in Fig. [2] It is well described by a Gaussian distribution 
with o = 0.09. That is, the model describes the mean behavior of SF as a function of 
luminosity, rest-frame wavelength and time to within ~10% of the measured values. 
However, the analysis based on two-epoch data cannot provide information about SF 
variance in a bin with fixed values of Mj, AIrf, and Xrf- This is because the SF is 
constructed using the magnitude difference measurements from all the quasars. Their 
scatter measures the mean value of the SF, but provides no information about the SF 
variance among individual objects. To measure the latter, individual light curves must 
be available. 

The distribution of SF t s /SF moi i e i values for the sample analyzed here, where SF is 
evaluated for every individual object from its light curve, is also shown in the upper left 
panel in Fig.[2l It is much wider than the corresponding distribution based on two-epoch 
sample (which covers a larger range in AtRp), clearly non-Gaussian, and well-fit by a sum 
of two Gaussians. In fact, the scatter in SF b s /SF m0C i e i for fixed values of Mi, AtRF, and 
Xrf, is similar to the scatter induced by the variance of these quantities. The published 
analysis, exemplified by VB, dVBWL and 104 work, captured the mean trends, but 
not the surprisingly large scatter around the mean behavior (o ~ 0.4). This new result 
demonstrates that there is significant additional information that can be extracted from 
well-sampled light curves for individual objects. 

At this stage of analysis it is not clear what causes the SF scatter. One obvious 
candidate is intrinsic stochasticity of the process causing variability (e.g. see discussion 
in [180 ) - We are also revisiting issues such as photometric calibration of SDSS stripe 82 
data and the robustness of SF analysis to non-Gaussian outliers (the distribution of two- 
epoch magnitude differences follows an exponential, rather than Gaussian, distribution; 
104). It will also be helpful in the future to repeat this analysis for better- sampled light 
curves. 
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FIGURE 2. Distribution of SF b s /SF mot i e i. The distribution from analysis by 104 based on two-epoch 
data is shown by the dashed curve (scaled by a factor of 0.2). The symbols with error bars show the 
distribution obtained here, that is based on SF measured for each individual object (n is defined as the 
number of points in a bin (Nbin) divided by the bin width (A&j n )> an d Njot is the total number of points 
used for each histogram; that is, the enclosed area is unity). Error bars are computed as \/Nbi n / (NroT^bin)- 
The curve passing through the data points is an empirical fit that is the sum of the two solid (Gaussian) 
curves below it, with relative normalizations of 44% and 56%. The same total distribution is shown as 
the dark histogram in the other three panels, where it is compared to analogous distributions for radio 
(top-right), X-ray (lower-left), and infrared (lower-right) subsamples. 



The SF for Subsamples with Detections at Other Wavelengths 

We investigated whether the distribution of SF t s /SF mo( iei for the sample analyzed 
here varies among various subsamples. First we compared the distribution for the whole 
sample to the distribution obtained for the apparently brightest 10% objects in the i band, 
and did not detect any statistically significant difference. We also created subsamples 
that are detected in 2MASS, FIRST, and ROSAT, as listed in [H. Out of the 10,370 total 
data points, 1856 (or ~400 quasars) have infrared (/, H, and K) detections, 766 (~150 
quasars) have radio detections, and 309 (~70 quasars) have X-ray detections. 

The distribution of SF £, s /SF moc t e i for each of these subsamples is compared to that of 
the entire sample in Fig. [2] While the ROSAT sample seems to follow the distribution for 
the whole sample, the FIRST and 2MASS samples appear skewed toward lower values 
of SF b s /SF moe j e [, indicating less variability compared to the optically selected sample. 
We also investigated the behavior of SF o t, s /SF m0c iei values as a function of optical and 
infrared colors / — K and J — K, optical-radio "color" i — t, and X-ray-optical "color" 



3.0 
2.5 
g2.0 

O 

bi 1 1.5 



C/l 



. ° 1.0 
0.5 
0.0 



3.0 
2.5 

g2.0 
o 

7j 1.5 
1.0 

/) 

0.5 
0.0 



, ■ ■ ■ 1 ■ , ■ ■ 1 ■ ^ ■ 1 ■ ■ ■ ■ 1 ■ ■ ■ ■ 1 r- 

sl'ope=0.00 ± 0.02, : 
\ . y-!nt = 0.82 ± 0.06 - 


—~ . . . 1 n • ~— 1 . . . . 1 . . . . 1 ~- 

• ' ' ■ slope=0.04 ± 0.03, : 
- ■ -'■ y-int = 0.77 ± 0.04 - 


• ■ ! • -J- ' ' • •* . 


yv'*$^^ y p ■ ■ - 


2.0 2.5 3.0 3.5 4.0 1.0 1.5 2.0 2.5 
I-K (mag) J-K (mag) 


slope=0.04 ± 0.01, : 
. ••• . y -int'=0.69 ± 0.05 : 


slope=0.03 ± 0.06, : 
: y-int=0.71 ±' 0.51 - 


I ■ . ■ I ■ ■ ■ I . ■ ■ I ■ ■ ■ I ■ ■ : 





4 6 
i-t (mag) 



6.5 



7.0 7.5 8.0 
x-l (mag) 



8.5 



FIGURE 3. The ratio of observed SF to model SF is plotted against i — K,J — K, i — t, and x — i colors as 
dots. These data points are then divided into bins of jV points, where is set to 10% of the total number of 
points. The median values for each bin along with 1 a error bars are shown, where a — 0.926(interquartile 
range) /\/N — 1. The over-plotted solid lines are linear regressions for the median values (the fit parameters 
are listed in the upper-right corners) and the over-plotted dashed lines are the same, but for the entire 
sample of 10,370 data points. 



x — i (t and x are radio and X-ray AB magnitudes, see [§(]). Figure [3] shows that the 
quantity SF b s /SF m0l j e i is independent of i — K and J — K colors. However, there seems 
to be a positive correlation with i — t: optical variability increases with radio loudness, 
in agreement with VB, who found that radio-bright quasars are about 1.3 times more 
variable. There is no significant correlation with x — i color. 
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